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Abstract 

Let p : N — > M be a surjective map of smooth manifolds. We are 
concerned with singular perturbation problems associated to a pair of 
second order positive definite differential operators with no zero order 
terms, that are intertwined by p. We discuss the associated random per- 
turbations of stochastic differential equations and present a number of 
examples including perturbation to geodesic flows and construction of a 
Brownian motion on S 2 through homogenisation of SDE's on the Hopf 
fibration. 
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1 Introduction 

The principal idea of singular perturbation is to deduce long term trends of 
a complex system from that of a relatively simple one for which some ob- 
servables are known. We are mainly concerned with random perturbation 
of stochastic differential equations with intertwined diffusion structures. A 
diffusion operator B on a smooth finite dimensional manifold A^ is a second 
order differential operator with positive definite symbol and vanishing zero 
order term. Its symbol a B is a real valued bilinear map on the cotangent 
bundle, determined by 

1 

a (dg h dg 2 ) = - {B{g 1 g 2 ) - {Bg x )g 2 - gi(Bg 2 )) , 
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where gi, g 2 are real-valued C 2 functions on N. Let n = dim(iV). In a local 
coordinate system, let 



where a it j and b k are smooth functions with the nxn matrix valued function 
(dij) positive symmetric. Then for y e N, 



The operator B is said to be elliptic if the symbol is strictly positive and to 
have constant rank if the rank of its symbol a® is constant in u. 

Let p : N — > M be a smooth onto map between smooth manifolds N 
and M. Denote by Tp : TN -» TM the differential of p. If B is an operator 
on N and A an operator on M such that 



for all real valued C 2 functions / on M, we say that B and A are intertwined 
or i3 is over A. The simplest intertwining is given by projection of a product 
space to one of the factor spaces. 

Let T u p : T U N — > T v(u )M be the differential of the map p at u, a linear 
map between tangent spaces. Its kernel is said to be a vertical tangent 
space, a subspace of the tangent space T U N and is denoted by VT U N. The 
vector bundle VTN = U u VT n N is called the vertical tangent bundle. We 
say that a diffusion operator B° on iV is vertical if B°(f o p) = for any C 2 
function / on M. 

Suppose that B is over A and that {a£, x e M} has constant rank then 
there is a unique smooth lifting map 



such that T u p o f) w is the identity map (Proposition 2.1.2 in Elworthy-LeJan- 
Li lfTTlO . The image t) M induces a smooth distribution, called the horizontal 
distribution associated to A. In the case of A elliptic let HT U N be the image 
of t) u then 



In the case where p : N — > M is a principal bundle this induces an Ehres- 
mann connection on iV with a corresponding connection 1-form. 





B(fop) = (Af)o P 



t) u : Image[^ n) ] C T p(u) M -> Image(0 C T M iV 
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Let A — \ YhLi LxiLx, + L Xo , where Xi are vector fields and L Xi de- 
notes Lie differentiation in the direction of X { . For short we also write 
A — \ YliLiX? + X . Suppose that X e Image o A we say A is cohesive. 
An elliptic operator is cohesive. The horizontal lift of A is: 

where Xi(u) is the horizontal lift of the tangent vector Xi(p(u)). By Theorem 
2.2.5 in mm, there is a unique vertical diffusion B v such that B = A H + B v . 
If u is a regular point of the map p, the B v diffusion starting from u stays 
in the sub-manifold The vertical diffusion operator is 'elliptic' if 

af : VT U N x VT U N ->■ R is strictly positive definite. 

Let A be a cohesive diffusion on M and B° a vertical diffusion. Let 
C e = \A H + B°, or C £ = A H + \B° where e is a real number which we 
take to zero. We would like to understand the asymptotic behaviour of the 
solutions of 4 = C e as e goes to zero. 

For simplicity write £ e = + £ x . Let us expand f e (t, y), a solution of 
the parabolic differential equation j^f e (t,y) = C e f e (t,y) in e: f e = ~/ + 
/i + e/2 + o(e). What can we say about f?. Is there a diffusion operator L 
such that fi solves = L? See e.g. the book of Arnold fl3]] in the context of 
perturbation to Hamiltonian systems and a recent book of Stuart-Paviliotis 
If26ll on multi-scale methods. 

We now introduce notations concerning Markov processes associated to 
diffusion operators. Let J 7 , Tt, P) be a filtered probability space. Let 
4>t(y, oj) be a family of strong Markov stochastic process with values in a 
manifold with 0o(y) = V- The ui variable will be suppressed for sim- 
plicity. Define a linear operator P t on the space of bounded functions by 
Ptfiy) = lo)), where E denotes integration with respect to P. Then 

(Pt, t > 0) is a semigroup of operators with P the identity map. Its infinites- 
imal generator C is defined by Cf = lim^o with domain the set of 
functions such that the limit exists. For such /, P t f solves the parabolic 
equation ^ = C with initial function /. On the other hand given any dif- 
fusion operator C, there is a strong Markov process whose infinitesimal 
generator is C. This can be seen by introducing a Hormander form repre- 
sentation of C and a stochastic differential equation. When C is reasonably 
smooth the Markov process is continuous in t for almost surely all uj G Vl 
and the Markov process is said to be a diffusion (process) with generator 
C. If C = \ A, where A is the Laplacian operator for a Riimannian metric, 
we say the diffusion process is a Brownian motion for that metric. 

The dynamic picture of the perturbation problem is as following. Let 
y\ be a diffusion operator with initial value y associated to C + eC\. The 
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y\ process follows roughly the orbit determined by £ with negligible dif- 
ferences. However on a large time scale of order 1/e, the deviation of the 
perturbed orbit from the unperturbed one becomes visible. In general we 
ask whether there is a function of y\ that does not change with time when 
e = 0. If so, denote this function by F. The variable F(yf) should vary 
slowly when e — > and in the limit F(y\) may converge to a Markov pro- 

cess whose probability distribution is determined by a diffusion operator L. 
It is desirable to find all conserved quantities and their explicit probability 
distribution, i.e. the limiting diffusion operator L, which is said to be the 
effective motion. There is extensive literature on this and we refer to the 
following books and the references therein: Bensoussan-Lions-Papanicolau 
and Freidlin-Wentzell HMD . 

We now consider the problem at the level of stochastic differential equa- 
tions. Every diffusion operator, if sufficiently smooth, can be represented 
as sum of squares of vector fields, the so called Hormander form repre- 
sentation. We write the two diffusion operators £ and Ci in Hormander 
form: C = \YZi x ? + x o, and A = § ££i Y? + Y . Let (6*,t$ be 
independent one dimensional Brownian motions on a given filtered proba- 
bility space J 7 , J^,P) with the usual assumptions. Let (jf t {y) denote the 
solution to the stochastic differential equations (SDE) driven by the vector 
fields -^Xi, \X Ql Yi, Y , i — 1, . . . , m, with initial point y: 

d Vt = 4^ Yl X *(yd ° db t + ~ e X o(y e t)dt + J2 Y M) ° dwi + Y (yl)dt. 

i j 

The solutions are continuous Markov processes (diffusion processes) with 
generator C e . 

Main Results. The structure of the paper is as following. In section 
2, we present several intertwining examples to illustrate the terminologies. 
In section [3] SDE's on the Hopf vibration are investigated. We construct 
stochastic process with generator a hypoelliptic horizontal diffusion on S 3 
from two vector fields one of which the vertical vector field induced by 
the circle action and the other a horizontal vector field induced from an 
element of the Lie algebra of norm 1 . In particular we construct a Brownian 
motion on 5" 2 through homogenisation. 

In section [4] the state space is the frame bundle of a Riemannian man- 
ifold. The frame bundle is closely related to the tangent bundle of the 
manifold, notably suitable first order differential equations on the frame 
bundle correspond to second order differential equations on the manifold. 
It is also the natural space to record the position and orientation of a par- 
ticle. For intertwined diffusion models on the frame bundle of a complete 
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Riemannian manifold M there are two basic slow motions: the horizontal 
diffusions and the vertical diffusions. The horizontal diffusions appear nat- 
urally in Malliavin calculus. The solution flow of a vertical diffusion can be 
considered as a random evolution of a linear frame and is an interesting 
object in geometry, e.g. it was examined in Brendle-Schoen QSj in connec- 
tion with the question as to whether positive isotropic curvature condition 
is preserved by R. Hamilton's ODE. 

In section 14.11 Perturbations to vertical diffusions on the orthonormal 
frame bundle, whose projection to the manifold is a fixed point are stud- 
ied. The effective motion on the manifold is a diffusion which is of the same 
'type' as the perturbation (Theorem 14. ID . Special care has to be taken in 
this case to avoid assumptions on the injectivity radius of the orthonor- 
mal frame bundle. In section 14.21 a suitable perturbation model to the 
geodesic flow, a first-order ODE on the frame bundle, is developed. The 
perturbation is of Ornstein-Uhlenbeck type and in the limit we see a dif- 
fusion with generator n(n 4 _ 1} A g , A H being the horizontal Laplacian, and a 
rescaled Brownian motion to the manifold. This study relates to approxi- 
mation of Brownian motions and by extension that of stochastic differential 
equations. We would also like to compare this to the philosophy in Bismut 
0, that x = h{—x + w) interpolates between classical Brownian motion 
(T — > 0) and the geodesic flow (T — > oo). In section [4741 perturbation to the 
the semi-elliptic horizontal flow is considered and we obtain an effective 
motion that is transversal to the holonomy bundle. 

In terms of ellipticity our SDEs have the following features. The un- 
perturbed system can be elliptic or hypoelliptic. The perturbations could 
be elliptic, hypoelliptic or degenerate. We used two types of scalings: the 
standard scaling and the scaling of 'Ornstein-Uhlenbeck type'. A related 
problem on commutation of linearisation with averaging is discussed in a 
paper in preparation II2T11 . 

2 Some Basic Examples 

Intertwined structures occur naturally. One such standard example is pro- 
jection of a product space. The other example is a principal bundle p : 
iV — >■ M with a group action G which acts freely. The latter introduces 
a twist in the product structure. These will include projections of groups 
to their quotient groups and that of manifolds to their moduli spaces. See 
[TTTTl and Liao [1221 for a discussion of diffusion processes on symmetric 
spaces. We give some examples which illustrate the procedure of averaging 
with intertwined structures and the local structure of frame bundles. Let 
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us begin with the trivial example of a cylinder K x S 1 . Let z denote the S 1 
direction, p(x, z) = z, A = J^, and B = sin z£ + J^. Then A H = J^. The 
effective motion associated to associated to hypoelliptic flows sin z-£ + 
converges to a simple completely degenerate motion. 

Example 1. Take iV = M 3 with the Heisenberg group structure. For 
(x, y, z) e M 3 define p(x, y, z) = (x, y). Let Y x = £ + \y-§- z and Y 2 = ^ - \x-§- z 
be the left invariant vector fields on iV associated to (1, 0, 0) and (0, 1, 0). 
Let A H = \{Y 2 + Y 2 2 ). Consider B° = \£> - f . Then C e = \{Y 2 + Y 2 2 ) + 
is over the cohesive operator A = \(-^ + Jp) on M 2 . The associated 
horizontal lift is (u,v) !->■ (u,v, \xv — \yu). The vector fields Y X ,Y 2 are 
respectively the horizontal lifts of ^ and -j^ and A H is the horizontal lift 
of A which is seen to have constant rank 2. The 'slow part' of the fully 
elliptic C e diffusion converges weakly to the A H diffusion, a hypo-elliptic 
diffusion, as expected. We illustrate what does it mean by the 'slow part' 
by the following trivial, but explicit, SDE: 

dx\ = cos(zl) o db\ — sin(^) o db 2 , dy\ = sin(^) o db\ + cos(^) o db 2 t 

dz t =\( x t sin(^) - y\ cos(4)) o db] + -(x e t cos(^) + y\ sin(^)) o db 2 

1 , 1 r , 
H — -pdw t z.dt. 

V e e 

Two obvious slow variables are {x\, yl). In our terminology the slow part is 
the solution to the following SDE parametrized by z\ : 

dx\ = cos(^) o db] — sin(^) o db 2 , dy\ = sm(z^) o db] + cos(^) o db 2 

1 1 
dz t =2^* sin (4) - ft c °s(4)) ° db t + g( £ t cos (4) + y| s in (4)) ° • 

The law of the process in independent of e and is \(-^2 + Jp-) — ygf^j + 

x af^ + l^ 2 + y 2 ^~&' Example 14. II for an example on the frame bundle 
that is in the same spirit. The third component of the Markov process 
associated to C e converges to the stochastic area of the limits of the first two 
components: | f* x e s dy e s — ~ J* yldx e s converges. This means taking stochastic 
area and taking e — > commute, as expected. 

This example can extend to the case of a general connection given by 
t)(x, y ,z)(u, v) = {u,v,r\u + r 2 v) allowing the functions to depend on z. 
Assuming that — ^) 2 is strictly positive, f) is determined by the oper- 
ator A H := \{£ + nf z ) 2 + + r 2 f z ) 2 and A, page 21 Elworthy-LeJan- 
LiUllI]. Let Xi = f- + r x -§- and X 2 = §- + r 2 f-. The first prolongation of 
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span{X!,X 2 }, i.e. span{XL,X 2 , [X 1; X 2 ]} has full rank at each point. Con- 
sider a function a such that the invariant measure jj, of (7 — r\ — r|)|- + a-j^ 
is the standard Gaussian measure. Let 7 be such that 7 — r\ — rf > c > 
some c. The + ^(7 — r\ — r|)|- + diffusion converges to an elliptic 
diffusion on the Heisenberg group if are not constants a.s. with respect 
to \i. 

Example 3. A non-relativistic quantum mechanical diffusion lives nat- 
urally in IR 3 x 5*0(3), the orthonormal frame bundle of IR 3 . Its spatial pro- 
jection lives in IR 3 . Studies associated to quantum mechanical equations, 
mainly the continuity equation describing the probability density of the 
quantum equation, have intertwined structures on p : R 3 x 50(3) — > IR 3 . I 
am grateful to D. Elworthy to bring my attention to the paper of Wallstrom 
11281 where limits of stochastic processes in IR 3 x 50(3) are discussed. The 
Bopp-Haag-equations have one free parameter I and its solutions converge 
to that of an equation with Pauli Hamiltonian as / — > 0. The Bopp-Haag- 
Dankel stochastic mechanical diffusions IR 3 x 5*0(3) were introduced by 
Dankel, describing a diffusion particle with definite position and orienta- 
tion. The Bopp-Haag -Dankel diffusions on IR 3 x 50(3) are given by a simple 
SDE with drift given by a Pauli spinor (solution of quantum equation asso- 
ciated with Pauli Hamiltonian with parameter /). In [28] it was shown that 
for spin \ wave functions and regular potentials the process parametrized 
by / converge to a Markovian process onto IR 3 , due to the averaging out 
of the orientational motion. The spatial projection describes the spatial 
motion of the particle without its orientation. 

Example 4. Let G = SO(ri) and tt : IR™ x G ->> IR n the projection to its 
first component. Let g = so(n) be the Lie algebra of the Lie group G. For 
each x e IR n , let h x : T x M. n ~ IR n — > g be a linear map varying smoothly in 
x. The map (x, v) i->- (x, h x (y)) can be considered as the horizontal lifting 
map through (x, I) where I is the identity matrix. This induces on IR n a 
non-trivial covariant differentiation V. Let e e IR n , consider the SDE 

dx t = e x g t o db t + eg t e dt 

p 

dg t = £ih x i(g t db t )g t + eh Xt (g t e)g t dt + V5^2z k (x t , g t ) ° dw\ + 5Z (x t , g t )dt. 

k=l 

where (b], w%) are independent 1 -dimensional Brownian motions, b t — (bj, . . . , 6"), 
w t = (wl, . . . , wl), and Z k : R n x G ->■ TG with Z k (x, g) G T g G. We denote by 
o Stratonovich integration, which must be used in the manifold setting and 
the correction term should be computed. When h = this corresponds to 
the flat connection. We consider three types of scalings: 1) 8 = 1, e\ = </e 
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and e — > 0; 2) ei = e = 1 and 5 — > 0. For the third type take ei = 0, e = 1 
and 5 — > oo. In case 1) it turns out that the solution x t is a slow variable, 
despite the involvement of g t . 

Let us now describe the model using the language of orthonormal frame 
bundles. Let v i-> t)(x, g )(v) be the horizontal lifting map. Then t)( X , g )(v) = 
t)(x,i)(v)g where / is the identity matrix. Denote by h x (v ) the lifting map at 
(x,/). Define 6( X , g )(v,w) = g' 1 ^). For (t>, A) e I" x g, define w (x j)(v,A) = 
A — h x {v) and 

m(x, g )(v, Ag) = g~ l W( xJ) {v, w)g = g^Ag - g~ l h x (y)g. 

The above example is the orthonormal frame bundle, R n x SO(n) — > R n , 
of R n with the group G = SO(n) acting on the right of R n x SO(n). This 
describes the local structure of the orthonormal frame bundle of a Rie- 
mannian manifold. The frame bundle of R n can also be represented as 
the Euclidean group. See section [4] for limiting theorems concerning the 
mentioned SDE in the context of a general manifold. 



3 Homogenisation on the Hopf Fibration 

Hopf fibration occurs in multiple situation in physics: in quantum systems 
and in mechanics, see e.g. Urbantke Il27ll for an account. Hopf fibration 
is the principle bundle tx : S 3 — > S 2 with S 1 acting on the right. Here 
S n denotes the n-sphere. As pointed out by M. Berger in 1962, this is a 
non-trivial collapsing manifold. The sphere S 3 is equipped with a metric 
inherited from R 4 . The collapsing was achieved by shrinking the length by 
a scale of e along the Hopf fibration direction and leaving the orthogonal 
directions unchanged. In a paper in preparation we study the dynamics 
associated to collapsing manifold Il20ll . 

It is convenient to consider the representation by unitary groups: S 3 
is identified with SU(2), S 1 with U(l) and S 2 with SU(2)/U(1). A typical 
element of SU{2) may be expressed as (z, w), where z, w G C are such that 

\z\^ + |w| 2 = 1, or as a matrix _ J . The right action by e %e e U(l) is 



w 

(z, w) i — y (e ie z, e tB w), which can be considered as right multiplication in the 
group SU(2) by elements of the form e 10 ~ (e td , 0): 



i 1 ^ i /in 

w z \ w z / \ 



z —w; \ ( z —w \ ( e ld 



-id 



The Hopf map tx : SU{2) — > S 2 is a submersion, 

n(z, w) = (Re(2zw), Im(2zw), \z\ 2 — \w\ 2 ). 
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The map T u tx can be better visualised if S 3 is treated as a subset of M 4 , 
writing z = y x + iy 2 , w = y 3 + iy 4 , y = (y u y 2 , y 3 , yd e R A , 

I 2/3 2/4, 2/1 2/2 \ 

Tyvr = 2 -y 4 y 3 y 2 , -yi 

\ 2/i, 2/2, -2/3, -2/4 / 

The vertical tangent spaces are the kernels of Ttx. It is easy to check that 
he vector field V(y 1 , y 2 , y 3 , yd = -j/ 2 #i + y\d 2 - y 4 d 3 + y 3 d 4 is vertical. Back 
to the principal bundle picture, V((z, w)) := (iz,iw) is the fundamental 
vertical vector field, associated to i in the Lie algebra of U(l). 

The Lie algebra su(2) is the set of matrices such that A + A T = and 
with zero trace: 

We take a real valued inner product on su(2), (A, B) := \ trace AB* , and the 
following orthonormal basis: 

Xl = ( o -i ) ' Xa = ( i o 1 ) ' Xs = ( i o ) • 

Note that X x is adjoint invariant under the circle action and so is the 
linear span of {X 2 ,X 3 }. Denote by X* the left invariant vector fields as- 
sociated to Xj. Let us define a distribution D = spanjX^Xg}, which is 
obviously left invariant with respect to the group action on S 3 . The span 
of the left invariant vector fields is also right invariant under the circle ac- 
tion. This is due to the fact that ue ie Xi e D uel e = u(Xie~ 2id )e ie e D u e id for 
i — 2, 3. Then T U S 3 = [kerT u 7r] © D u defines an Ehresmann connection on 
the principal bundle and a horizontal lifting map. 

Let V L be the left invariant linear connection and V the Levi-Civita con- 
nection for the bi-invariant Riemannian metric on the Lie group SU{2). De- 
note by A the Laplacian on S 3 . Let Ah = Ei= 2 V L c/f (X*, X*) = ^ =2 C Xi L Xi 
be the hypoelliptic Laplacian corresponding to the Horizontal distribution 
generated by the left invariant vector fields {X 2 , Xg}. 

3.1 Construction of Brownian Motion on S 2 by Homogenisation 

Let Y to be vector in span{X 2 ,X 3 }. Since Y e tB remains in span{X 2 ,X 3 } 
the SDE du\ = u € t Y gldt + -^al\X x o db t on S 3 makes sense for any stochastic 
process g\ e U(l). Note that {X 1 ,X 2 ,X 3 } is a Milnor frame flZl with 
structural constants (—2, —2, —2), 

[X 1; X 2 ] = -2X 3 , [X 2 , X 3 ] = -2X 1; [X 3 , Xi] = -2X 2 . 
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UY = c 2 X 2 + c 3 X 3 ^ 0, span{X 1 ,Y ,[Y ,X 1 ]} = span^, X 2 , X 3 }. By 
the structural equations [F ,^i] = 2c 2 X 3 - 2c 3 X 2 and {X 1 ,Y , [l^Xx]} is 
linearly independent following from the non-degeneracy of the matrix 

1 

c 2 -2c 3 
c 3 2c 2 

It follows that the SDE under discussion is hypoelliptic. 

The Hopf map n projects a curve u t in SU(2) to one in S 2 . A curve x t in 
S 2 lifts to a horizontal curve x t in SU (2) through the horizontal lifting map 
induced by the Ehresmann connection. 

Theorem 3.1 Let (b t ) be a one dimensional Brownian motion and take u e 
SU(2). Let (u e t , gl) be the solution to the following SDE on SU(2) x U(l), with 
u e = u and g e = I, 

1 1 

du\ = u e t Y gldt H — -=.u\Xx o db t , dg\ = —=g\X\ o db t . 

V e V e 

Let = 7t(m^) and x ^ its horizontal lift. Then x\ converges in probability 

to the hypoelliptic diffusion with generator CF = ^\Y \ 2 A H . IfY is a unit 
vector, x\ converges in law to the Brownian motion on S 2 . 

Proof Let a\ e S 1 be such that u\ = x\a\ where x\ is the horizontal lift of 
x\ through u using the connection determined by {X£,X£}. Then a e = 1 
and 

d~x\ = TR (aV -i o du\ + ((air 1 o datT(ul). 

All stochastic integration involved in the above equation are Stratonovich 
integrals. Thus 



dx\ = TR Kr i [u e t Y g$dt + ~jf t X x o db t j + ((a\) o d^T 1 )*^)- 

Since x\ is horizontal, u(dx e t ) = 0, we obtain 

(a e t ) o d(a\Y 1 = -w i% (^ulX^aiy 1 o db^j = --^=o|Xi(<$ _1 o db t . 

It follows that d(a e t y 1 = -^a^i(a t e ) _1 o db t , (a\) = (gl), and 

dx\ = u e t Y dt + -^=uJXi($) _1 o d& t - -^=xJaJXi(aJ) -1 o c#> t = x £ t g £ t Y dt. 
V e V e 



HOMOGENISATION ON THE HOPF FlBRATION 



11 



Since there is no Stratonovich correction term for dg t = g t X 1 o db t , the 
corresponding infinitesimal generator is |A 5 i where A 5 i is the Laplacian 
on S 1 . Let F : S 3 -» R be any smooth function. Since Y e span{X 2 , X 3 }, 

F(x e t ) = F(u )+ dF(x e s glY )ds = F(u ) + y2 dF{xlX j ){x e s X j x € s g^)ds 
Jo j =2 Jo 

= F(u ) + J2 f dFix^X^Xj, g*Y )ds. 

3=2 J ° 

The two real valued functions on S 1 , g i-> (X 2 , gY ) and g i-> (X 3 , gY ), are 
eigenfunctions of Agi. Then for j = 2, 3, 

dFixlXj^Xj, g$Y ) - dF{u G Xj){X h Y ) 

= V f v L rfF(^x fc , %-)<x is ^y ) (x fe , ^y )ds + 4= /* dF^x^Xj, glx^dK 

k=2 ^° ^ 6 ^° 

+ -/ ^^(x^^yo)^. 

e </ 

Applying the identities X? = —I we obtain 

3 

= F(tio) - e ^ (dFixtXjXXj, g e t Y ) - dFiuoX^Xj, Y )) 

3=2 

+ e Yl f^ L dF(xlX k) ff s X J )(X v glY Q ){X k) glY Q )ds (3.1) 
+ V^Y,[ dF(x e s X j )(X j ,gtX 1 Y )db s . 

3=2 J ° 

Since F is a smooth function on compact manifolds, the probability dis- 
tribution of {x\,e > 0} is tight, see Lemma [3721 below. We now move to 

the canonical probability space with the standard filtration F t . By Lemma 
13.31 below, conditioning on the filtration F s , the canonical filtration on the 
canonical probability space, 

3 4 

e E r v L c /F(^x fc ,^x J )(x„^y )(x fc ,^y )^ 

j,k=2 J l 

converges to 

V / V^FiuX^uXj) [ (X j ,gY )(X k ,gY )dgds. 

j,k=2 Jr J s" 
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Here dg is the Haar measure on S 1 . It is easy to check that 

/ (X 2 ,gY )(X 3 ,gY )dg = 0, 
Js 1 

either by direct computation or note that there is g' e U(l) such that 
g'X 2 = —X 2 and g'X 3 = X 3 and using the translation invariance of the 
Haar measure. Since there is an element of S 1 that maps X 2 to X 3 , 

[ (X 2 ,gY ) 2 dg= [ (X 3 ,gY ) 2 dg. 
Js 1 Js 1 

Note that 




2 



We conclude that x\ converges in distribution and its law is determined 
by the generator CF(u) = \\Y \ 2 V L dF(uX 2 ,uX 2 ) + \\Y Q \ 2 V L dF(uX 3 , uX 3 ). 
Since V L X* = for i = 2,3, £j =2 V L d(f o n)(Y*,Y*) = trace Vdf. Note 
also the Riemannian metric on S 2 is that induced from S* 3 , the process x\ 
has generator ||Y | 2 As2 and is a Brownian motion when Y is a unit vector. 

□ 



Lemma 3.2 Let fi e be the probability distributions of the stochastic processes 
(x\,t > 0) from the theorem. Then {//, e > 0} is relatively compact. 

Proof Write y\ = x\ for simplicity. Let [i n be a subsequence from 

corresponding to a sequence of numbers e n . We wish to prove that it has a 
weakly convergent subsequence. It is sufficient to prove that the family of 
measures /i n is tight, i.e. for every 5 > there exists a compact set Ks C M 
such that n n {Ks) > 1 — 5 for all n. As probability measures on the space 
of continuous paths on M, /i n (a : a(0) = y ) — 1 where a : [0, 1] ->■ M is 
a continuous path on M. For any yi,?/ 2 e M, let : \I x \l > } ; i be a 
smooth function that agrees with the Riemannian distance function when 
d(yi, y 2 ) < a/2 where a is the injectivity radius of M and <f>(yi, y 2 ) = 1 when 
d(yi,y 2 ) > 2a. This is possible by taking = a o d where a : M + — >■ R is a 
suitable bump function with a the identity function on [0, a/2]. Then is a 
distance function on M that generates the same topology as d. The family 
of measures {/!„} is tight if for any a, i] > there exists < 5 < 1 such that 
there is an e > 0, with 

P w: sup <p(y e s n , y e t n ) > a) <i], when e < e . 

V \s-t\<6 I 
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In the proof of the Theorem, take F(y) = <p(y, u). Then by formula (I3.1D . 

Esup 2 (^, u) < (p 2 (y e , u) + Ce + e5 

s<5 

for come constant C. Let (j) e t (y,uj)) denote y e (u) with y e Q {u) = y. Let 9 S 
denotes the shift operator in the Wiener space. By the Cocycle property, for 

s < t, 

E sup (j> 2 (y e s \yt n ) = EE{ sup <f>\z, cf> £ t _ s (z, 6 t - s (cu))M n = z} < Ce + e6 

\s-t\<5 \s-t\<6 

and the required tightness holds. □ 

Let y™ be a family of Markov processes on a Riemannian manifold M 
that is relatively compact. Represent this as the coordinate process on path 
space with measure fi n , the distribution of y e . Suppose that \i n converges 
weakly to fi. Let F : M — > R be a smooth function with compact support. 
Let A be a diffusion operator. Suppose that 



F(Xi) - /■'( A'n) - / AF(X r )dr W///„ -> 



for any function / that is measurable with respect to T a where (J 7 ,, s > 0) 
is the canonical filtration. Then Ji is the probability distribution of a A- 
diffusion. In fact letting Mf = F(X t ) - F(X ) - f* AF(X r )dr, then Mf is 
a jj martingale and /i is the solution to the martingale problem associated 
to A. The following lemma reflects this philosophy, let e n be a sequence 
converges to zero. We are interested in the term 

3 * 

e E / £ v L rfF(x:x fc ,x:x J )(x,,^y )(x fc ,^y ) c / s . 

Lemma 3.3 Let {y e t , hf) be a family of SU(2) x U (1) valued stochastic processes 
on a probability space such that the law of(y^, hf) agrees with that of{x\,gl) 

in Theorem \3. 1 1 Let (x\, g e t n ) be a weakly convergent sequence. Let := y\ n . 

We may assume that y n converges to y. almost surely. Let F : SU (2) — > R be 
a smooth function. Define 

CF(u)= y2v L dF(uX kl uX 3 ) [ (X 3 ,gY )(X k ,gY )dg, 

3 

A e F(yt,gl) = V L dF(yl e X kl yl e X )(X v glY )(X kl glY ). 

j,k=2 
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Then the following convergence holds in L 1 , 

e A e F(y^g*)dr-+ J CF(y e r )dr, 
and for any real valued bounded function 4> on the path space, 

mVr,r < s) (F(y?) - F(y r :) - jf CFiy^dr^j -> 0. 
Proof By formula (15TTD . 

■ / - i=2 j=2 

+ e j* A n F{yl,gl)dr- J CF{y n r )dr. 

It is sufficient to prove that 

e[ e A n F(y e r ,gt)dr^ J CF{y")dr. 

Let t = s < ti < • ■ ■ < t n = t be a division of [s, t] with appropriate scale. 
Let Ai t = £ i+1 — £j. Assume that At.; = ^/e so ^ = ^= is large. On each 

interval ^-1, 



|V L rfF(^ e X fc ,^ e X,) - V L dF(^X fc ,^ e X,)| < C|^ e - y e t J ~ o^At,), 

where ~ means in the order of, after taking expectations. Let dg be the 
Haar measure on S 1 . It is the invariant measure for g e tt where g\ is solution 
to dg\ = ~g\X\dt. By Birkhoff's ergodic theorem on S 1 we obtain, 

t 

e £ V L rfF(^ e X fc , ylXJiXj, gfY ) (X k , glY )ds 

e 

~ £ Ati\7 L dF(yl i€ Xk, y^Xj) f (X j>9 Y ) (X k , gY )dg 

J s 



' V L dF(f re X k , y e re Xj)dr / (X,-, gY ) (X k , gY )dg 

Js 1 

V L dF(y;X kj y;X 3 )dr [ (X J} gY )(X k , gY )dg. 



□ 
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4 Perturbed Systems On Principal Bundles 

Let M be a smooth finite dimensional manifold. The frame bundle tx : 
FM — v M is a principal bundle with group action GL(n, M). Its total space 
is the collection of all linear isomorphisms u : IR n — >• T^M. Given a 
Riemannian metric on M, the orthonormal frame bundle 7r : OM — > M 
is a reduced bundle with group action 0(n) and the fibre at u : IR n — >■ 
T n ^M consisting of isometric linear maps. The total space of the frame 
bundle or the orthonormal frame bundle is a manifold in its own right. If 
M is oriented the orthonormal frame bundle consists of two components 
in which case we only need to consider the action by the component SO{n) 
of the group that contains the identity. We assume that n > 1. In all cases 
the group will be denoted by G, its Lie algebra by g and the right action by 
a e G is denoted by R a . In case of G = SO(n), q = so(n) is the space of 
skew symmetric matrices. 

Denote by TOM the tangent space of OM and by VT u OM the naturally 
defined vertical sub-bundle, VT u OM = ker[T u 7r]. If A belongs to the Lie 
algebra so(n), denote by A* the fundamental vertical vector field on OM 
induced by right multiplication, 

A *( u "> = ^i\t=ouexp(tA). 
at 

Then an o.n.b of so(n) induces a family of vertical vector fields that spans 
VTOM and VTOM is an integrable sub-bundle of the tangent bundle . 

A connection V on the tangent space of M induces a splitting of the 
tangent spaces of T u OM: 

T u OM = iri„()M © VT u OM. 

Let HTOM = U u HT u OM and VTOM = U u VT u OM. Then HTOM is a 
right invariant distribution and the splitting is an Ehresmann connection of 
TOM. Tangent vectors or vector fields are called horizontal (respectively 
vertical) if they take vales in TOM (respectively in VTOM). This deter- 
mines a linear connection V on M, a connection 1-form zu u : T u OM — > 
so(n), and a horizontal lifting map fj u : T n(u) M — > T u OM. 

To each e e M™, there is associated a standard (or basic) horizontal 
vector field H(e) given by u (->■ H{u){e) = f) u (ue). For A e so(n), [H(e), A*] 
is a horizontal vector field. For e,e e R n the vertical part of [H(e), H(e)] is 
given by the curvature form fl, [H(e), H(e)] = 2Q(H(e), H(e)). Let {e^} be 
an orthonormal basis of R™ and define Hi = H(u)(ei). 

Let (Q, J 7 , T t , P) be a filtered probability space with the usual assump- 
tions and let {w{, b l t , 1 < j < p, 1 < I < m} be independent one dimensional 
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Brownian motions. Let w t = (w\, . . . , w™) and b t — (b\, . . . , Let {X z , I = 
0, 1, 2, . . . , m} be a family of horizontal vector fields and {Zj, j — 0, 1, . . .p} 
a family of vertical vector fields. Consider the SDE 

f du\ = yfe x YZi X K M D ° dh \ + eX (uf)dt + V6 £J =1 Z 3 (u\) o dw{ + 8Z (u e t )dt, 
\ u € = u . 

(4.1) 

here e±,e and 5 are parameters. The infinitesimal generator of the SDE is 

h Eh X * + eX + 5 £™ =1 ^ + <^c 

4. 1 Perturbation to Vertical Flows 

A Riemannian connection V on M is a connection that is compatible with 
the Riemannian metric, with possibly a non-vanishing torsion T. We take 
the horizontal bundle on the principle bundle induced by this connection. 
We will assume that the connection is complete, i.e. every geodesic extends 
to all finite time parameter. This is so if every standard horizontal vector 
field is complete. Let w : T u OM — > so(n) be the connection 1-form, cor- 
responding to the given Riemannian connection V, which is determined 
by adjoint invariance and its values on fundamental vertical vector fields: 
(R g )*w = adig-^zu and tu(A*) = A. Let 6 U : T u OM W 1 be the canonical 
1-form such that 6 u (fy u (ue)) = e. Let V be the direct sum connection on 
TOM. For any vector v e T u OM and vector field U on OM, 

V V U = w- l d{w(U))(v) + 6- l d(0(U)){v). 

This connection V has zero curvature and a non-vanishing torsion in gen- 
eral. A formula for the torsion of V is given in Li lf2~TTl . 
In (14.11) take e x = e and 6=1. Let If = £ + ed where 

^ m \ P 

1=1 0=1 
Theorem 4.1 Assume that M has positive injectivity radius, {tu u [Zj(u)]}Y = i 
spans q, and the vector fields {X h l > 0} and {|V Xi X ; |,/ > 1} have linear 
growth. Let u\ he a solution with initial value u G OM, to the SDE 

m p 

du e t = ^Y^ ° db t + eX (u e t )dt + z M t ) o dw{Z (u e t )dt. (4.2) 

Let x\ = ir(u e t ) and x\ its horizontal lift. Let [i u he the invariant measure of 
the following SDE on G: 

in 

dg t = 2_ j TL gt w[Z j (ugt)'\ o dw{ + TL gt w[Z (ug t )]dt. 
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b(u) = J ^ VxjXj(u^) + Xo(ug) j d^„(£) 
aij(«) = / S ^{TR' g l ^i{ug),H, i {u)){TR' g l ^i{ug),H j {u)) d^iuig), 



g i=i 



Then x\ converges weakly with limiting generator C. For F : OM — >■ 
smooth with compact support, 



1 P 

CF(u) = dF(b(u)) + - a itj (uyVdF(Hi(u), H^u)). 



i,3=l 

Proof Since 5^ and belong to the same fibre we may define g\ e G by 

u\ = x\g\. If a t is a C 1 curve in G 

-r;\tua t = — ua t a t a r+t = (a t a t ) (ua t ). 

at dr\ r=0 

It follows that 

du\ = TR gl d~x\ + (TL (gt) -idg$y(uf). 

Since right translation of horizontal vectors are horizontal, zu(dul) = TL^-idgl 
and 

m 

dg\ = TL gt w[Zj(xlg e t )] o dw{ + TL g% w\Z^x\g^)\dt. (4.3) 

3=1 

By Ito's formula, dx\ = v^ELi Tvr(X/(M t e )) o dh\ + eTn(Xo(u e t ))dt so 

v 

d~x\ = t)* t (odxl) = V^J2 ^[TniX^))] o db\ + efjx«[T7r(Xo(« t e ))]d*. 

By assumption on the vector fields X;, the above SDE does not explode and 
n(ul) exists for all time. In terms of the group action, we have 

v 

dx\ = V^J2 M^f^O^t)] ° db\ + e^|[T i|9| 7r(Xo(x t e ^)]d*. (4.4) 
z=i 

Let /x e be the laws of the x^. We first show that {/i e } is tight. By Prohorov's 
theorem a family of probability measures is tight if it is relatively compact. 
Since x^ = u it suffices to estimate the modulus of continuity and show 
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that for all positive numbers a, rj, there exists 5 > such that for all e 
reasonably small, see Billingsley []6l] Ethier- Kurtz lfT3H . 

P(u : sup d(xl,x e s ) > a) < r). 

\s-t\<8 

Here d denotes a distance function on OM. The Riemannian distance func- 
tion is not smooth on the cut locus. The cut locus of OM is in general not 
predictable by that of M. To avoid any assumption on the cut locus of OM 
we construct a new distance function that preserves the topology of OM. 

Let x G M and 2a the minimum of 1 and the injectivity radius of M. 
Let cf) : R + — > M + be a smooth concave function such that 0(r) = r when 
r < a and 0(r) = 1 when r > 2a, e.g. 4> is the convolution of min(l, r) with 
a standard mollifier supported in the set {r : \r — y| < a/2}. Let p and p 
be respectively the Riemannian distance on M and OM. Then (pop and 
d := 4>o p are distance functions. For m g 7r _1 (x), 

o p( M , ~x\) ={<p o p)( M , ^) + f d(<P o p) ( V fe. [Ttt( W s ))] o cLB< 

■ / ° V 1=1 

+ / ed(0op)h-[T7r(Xo«))]d S 



o 



o p)( M , 4) + J Q d((l>o p) i^fe ^[Ttt(XK^))]^ 
/ Vd(0op)(T7r(XK<)),T7r(XK<))) ds 



i=i 



+ e jf d(0 o p) f 1 V T ^)(Tvr o X0«) + Tvr(Xo( M :)) J ds. 

Since o p has compact support and the vector fields concerned have linear 
growth, |Ttt(X,(i#)| < C(l + < [C + Cp~(x e s ,u e s )] + Cp(u,x e s ) some 

u G OM. The quantity C + Cp(x e s , u e s ) is bounded from the compactness of 
G and it follows that E[0 o p(w, x\)f) < Ci(t)((0 o p) 2 («, x%) + et) for some 
constant C. By the Markov property and the estimates below the required 
tightness follows, 

E[0op(51,^ t )] 2 ) < C x \t-s\. 
By the right invariance of the horizontal lift, 



^[T £Ss |7r(Xj(x^)] = TR^ej-iXiiul). 
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Let F : OM — > K be a smooth function. For V, the canonical direct sum 
connection on OM associated to V, 



p ft 

F(x\) = F(u ) + V~eJ2 dF ( TR ( 9 tr^i«)) dB l 

1=1 Jo 

+ 2 e S / ^dF^^yiX^ullTR^-iXtiu^ds 
1=1 ^° 

+ 2 e Z) / ^(Vx,X,(«p + Xo(«S)j£fa. 

By tightness and Prohorov's theorem we may take a sequence e n — >■ with 
the property that xt 1 converges in law to a probability measure fi. Let X 

be the coordinate process on the path space and Q t = a{X s : < s < t}. 
Since G is compact the following term has at most quadratic growth, 

(TR g -iXi(ug), Hi(u))(TR g -iXi(ug), Hjiu^^dg), 



and by the same argument J G Ym=i VXi(X[)(ug)[i u (dg) has linear growth. 
To identify the limiting process it suffices to show that for all real-valued 
smooth functions F on OM with compact support, 



CF(X r )dr I g d{i 



F(X t ) - F(X S ) 



converges to zero where g is any real-valued bounded C? s -measurable func- 
tion on the Wiener space and X t the canonical process. 

Let z™ be a sequence of random variables whose law agrees with that 
of x e j_ for some sequence e n and converges almost surely. Let g be a 

{z™, s < t}-adapted bounded function. For t > s, 



Eg F(z?) - F{z n s ) 



lF{z n r )dr 



E 



g / (A 6 "F - CF)(z?)dr 



->0, 



where A £n F is given by the bounded variation part in the formula for F(x e t ). 
The convergence holds since G is compact and also the invariant measure 
/i G for the elliptic SDE (14.31) is ergodic. The proof is standard and follows 
from the Lemma below. See e.g. Hasminskii lfT5l . Papanicolaou- Stroock- 
Varadhan Q. □ 



Perturbed Systems On Principal Bundles 



20 



Lemma 4.2 Let f be a bounded function with bounded derivative then 



The proof is completely analogous of that of Lemma 3.2 in II19II . In sub- 
intervals whose length is very small compared to 1/e we consider x e s as 
constants, and apply the ergodic theorem on each interval. With the size 
of the sub-intervals chosen correctly, the sum over all sub-intervals of the 
limits forms a Riemann sum. The convergence follows from the Cocycle 
property of the flows, estimates for the rate of convergence in the ergodic 
theorem and the regularity of the function A e f. 

In the theorem above the assumption on the injectivity radius can be re- 
moved in the case of the projection being a Brownian motion with bounded 
drift. See e.g. the estimates in [fLSll . We look into two special cases, when 
the horizontal vector fields are either right invariant (lifts of vector fields 
on the manifold M) or the standard horizontal vector fields. 

Example 4.1 (The Right Invariant Case) Let X h I = 0, 1, 2, ... m, be vec- 
tor fields on M. Define Xi(u) = t) u (Xi(n(u))) and we have 



The projection ir(u e t ) satisfies dx\ = \fl Ym=i Xi(x\) o db\ + eX (x e t )dt. For 
all e, x\ are \ V L Xi L Xi + L Xo -diffusions. The horizontal lift xt of xt are 



\ L^L^ + L Xo -diffusions. 

Example 4.2 (The Rotational Invariant Case) Let {q}" =1 be an o.n.b. of 

M", e e W 1 . Let H t (u) = H(u)(e t ), and H (u) = H(u)(e ) be horizontal 
vector fields. We have 




where (Esup s < t \R(f, e, s, t)f) " < C(t)ez for any > 1. 



p 



du £ t = ^Y^ Mu e t ) o db\ + eX (u € t )dt + Zj(u e t ) o dw{ + Z (u £ t )dt. 
i=i j=i 



n 



in 



dul 



y/e ^2 H i( u t) ° db t + eH (ut)dt + Z^u\) o dw{ + Z (u e t )dt. 

1=1 j=l 



Write xt 



l . Then 



dx\ = yfeH(xl)(jgl o db t ) + eH(x e t )(gte )dt. 



(4.5) 
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Its 'formal' Stratonovitch correction term vanishes. If x e = u then g\ is the 
identity matrix. Write dw\ = dt and let Zj = a{A* where {A,} is an o.n.b 
of so(n). Then 

dgt = ^oi(uf)gtAjodwt. 

If a 3 h are constants the process g\ is independent of e and x\ is a Markov 
process on OM. If furthermore e = 0, the law of x\, and hence that 

of x|, is independent of e. This follows from the independence of g\ and 

Finally x\ is a horizontal Brownian motion with projection x t a Markov 
process and a Brownian motion on M. This is the construction of Brownian 
motions of Eells-Elworthy [10]. The invariance is no longer true for e ^ 0. 
Remark: More generally if {& t (u)} is a family of Markov processes on 

OM with the property that ®t(ug) = &t(u)ipt(g) for some ip t (g) G G and 
a{7r($ r (M))|r < s} = cr{$ r (/u) : r < s}, then 7r($t(w)) is a Markov process. 
Denote by Q t (u ,du) the law of $t(«o) and let / : M —> R be a Borel 
measurable function, x t = 7r($ t (tt )), 

E{/(x t )|cr{x r , r < s}} = j(J o 7r)(«)Q t _ s (x s , du). 

It follows that f(foTr)(u)Q t _ s (x s ,du)= f(foir)(uil) s (g))Q t - s (x s g,du) = f(fo 
7r)(u)Q t -s(x s g,du). So E{/(x t )|cr{a; 7 ., r < s}} depends only on x s = ir(x s ). 
When e = 0, the flow of (I4.5D satisfies the rotational invariance condition 
and the horizontal lift of x t is a function of the path (x r , r < t). 

Example 4.3 Let a : M x I" 4 K" be a smooth map so that a(x) G 

L(R n ; R n ). Let{ ei }^ =1 be an o.n.b. of R n , e G R". Consider 

n m 

du\ = ^^t) u [a(n(u))ei]odb\+et) u [a(%(u))e ](ul)dt+'^2 Zj(u e t )odw J t +Z (u e t )dt. 

1=1 3 = 1 

The projection x\ = n(u e t ) satisfies: 

n 

dx\ = y/e ula(xl)(ei)odb l t +eula(x t j.)(eo)dt = \/eula(xl)odb t +eula(x < j.)(e )dt. 
i=i 

Let x\ be the horizontal lifting map of x\ and g\ be an element of G deter- 
mined by u\ = x\g\. Then dx\ = ^/eH(xf)g^a(xl) o db t + eH(x € t )g1a(x € t )(e )dt. 
When a(x) is not trivial the bounded variation term for f(x t ), where / : 
M — > R is a smooth function, will involve Yl% Vdf(ula(xf)ei, w|a(x|)ej) 
which is no longer a trace. It will also involve the derivative of X/. In 
this case it is useful to consider the system as perturbation of the vertical 
SDE about which we know a lot more. 
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4.2 Perturbation of Ornstein-Uhlenbeck "type 

We now describe the relation between horizontal equations on frame bun- 
dles and geodesic flows. Let P = GL(M) be the linear frame bundle over 
M. A vector field on a frame bundle can be considered as a second order 
differential equation on the underlying manifold as below. Fix e G R™ and 
H the isotropy group at e of the action G = GL(n, R) on R™. The tangent 
bundle TM can be considered as a fibre bundle associated with the princi- 
pal fibre bundle P with fibre R" . The total space E is P x R n / ~ where the 
equivalent class is determined by [w,e] ~ [ug~ x ,ge\, any g G G. Elements 
of the form where g £ H belong to the same equivalence class. It can 
be identified with the quotient bundle Pj H, whose element containing u is 
the equivalence class of the form {ug,g G H}. Denote by £ the coset H. 
Let a be the associated map: 

a eo : u G P — > ueo G TM. 

This induces a map w G T U P — > T u a eo (w) G T ueo TM. Each element v G TM 
has a representation d = ue , where u is unique up to right translation by 
elements of H. Furthermore a right invariant vector field W on P induces 
a vector field on TM. In fact if v = ue' = ue there is g G G with u' = ug 

and e' = flf -1 ^. Since a eo (w) = a e > Q {R g u), 

T u a eo (W(u)) = T u ,a e[ TR g (W(u)) = T u ,a e > o W(u'). 

This map W G TTP X w G TTTM is independent of the choices of e . 
Fix e . Any vector field that is invariant by right translations of elements 
of H induces a vector field on TM. Consider a horizontal distribution 
determined by a connection on TM and let W(u) = H u (e ) be the funda- 
mental horizontal vector field associated to e , the induced vector field is a 
geodesic spray X, i.e. in local co-ordinates X(x,v) = (x, v, v, Z(x, v)) and 
Z(x,sv) = s 2 Z(s,v), which corresponds to the geodesic flow equation on 
TM: 

dv t = -T at (vt)(v t ), & t = v t , cr(0) = 7r(u), v(0) = ue . 

Here T denotes the Christoffel symbol. The corresponding horizontal flow 
on P is given by u t = H(u t )(e Q ). 

Based on E. Nelson's Ornstein-Uhlenbeck theory of Brownian motions 
Il24l we ask the following question. What happens if we replace the driv- 
ing Brownian motion dw t by v t dt where v t is an Ornstein-Uhlenbeck pro- 
cess? Consider the position process z t in R™ with velocity process satisfy 
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the Langevin equation: 



dvf = — vl dt H — dw t 
e e 



where w t is a Brownian motion with values in W 1 and z = 0. The z| 
process converges to w 4 as e — >■ 0. The convergence holds almost surely and 
in fact the result holds if w t is replaced by any continuous function. We now 
interpret the convergence in terms of homogenisation. First we rescale the 
variables in space and time and setting v t = \Jev t , z t = \fez t . It is easy to 
see that z\ is the slow variable and i/ezt converges to a Brownian motion. 
In fact 

dv\ = — v\dt H — i=dw t , z\ = v\. 
e V e 

We must take care to take this model to the orthonormal bundle. First 
we are not allowed to rescale variables in non-linear spaces. We should 
not rescale the frame variable, in the orthonormal frame bundle, in space 
either. 

We shall consider the orthonormal frame bundle of an n-dimensional 
manifold with group action by G = O(n) or G = SO(n) if the manifold 
is oriented so that OM is a connected manifold of its own right. In the 
latter case we assume that n > 1. Let e G M n be a unit vector and 
{A k ,k = 1,2, ...,N = n{n — l)/2} be elements of g, and A G q. Let 
A* k be the corresponding fundamental vertical vector field corresponding to 
A k . Consider 

1 - 1 

du\ = H(u £ t )(e )dt + —= V A\{u\) o dw k t + -A* (u e t )dt. 

v e k =i e 

For 'e = oo', the equation can be considered as the 'geodesic flow' equation, 
as explained earlier. If x\ = n(u\) then x\ = u e t e . Note that the change of 
the velocity of the motion on M is always unitary. Due to the fast rotation, 
the geodesic has rapid changing directions and we expect to see a jittering 
motion and indeed we obtain a scaled Brownian motion in the limit if the 
rotational motion is elliptic. 

A related theorem is given in Dowell [J9l] stating that an Ornstein-Uhlenbeck 
position process on 2-uniformly smooth Banach manifolds converges. Those 
are manifolds modelled on 2-uniformly smooth Banach spaces. By a 2- 
uniformly smooth Banach space B we mean one with the property that 
there is a constant C > such that \ \x + y\\ 2 + \ \x — y\\ 2 < 2||x|| 2 + C||?/|| 2 
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holds for all x, y G B. The iterated Ornstein-Uhlenbeck processes in []9]] 
and the settings of manifolds with Lorenzo metrics, see e.g. Bailleul Ipfll, 
are also worth exploring. We expect interesting results arise for processes 
with infinite-dimensional noise. For a related work, central limit theorem 
for geodesic flows, we refer to Enriquez-Franchi-LeJan [121 . 

Let A* k denote also the left invariant vector field on G induced by A k G g: 
A* k (g) = gA k . Let A = \ J2 k ( A l) 2 + A o- Denote by A L the left invariant 
Laplacian on G. For all e, let u e = u and x = 7r(w ). Let {e^} be an 
orthonormal basis of R n and we may take e x — e . 

Theorem 4.3 Let M be a compact Riemannian manifold, u G OM. Let u\ 
be the solution to the SDE on OM: 

1 N 1 
du\ = H(u e t )(e )dt + — V o + -^*(w t e )dt, = u . (4.6) 

v e fe= i e 

//^4 is elliptic then ir(u\) and its horizontal lift converges in law. If further- 
more A = |A L then n(u\) converges in law to a rescaled Brownian motion 
with generator ^^^ A. Its horizontal lift converges in law to the diffusion 
process on OM with generator n ^_^ A H . 

Proof Define the Lie group valued process g\ by u\ = x\g\ and g e Q = I, the 
unit matrix. Following earlier computations and using the fact that A* is 
right invariant and {R g )*w = ad^ 1 )^, 



As in the proof of the previous theorem, we see that the family {x\} is tight 

and that it converges in law as e — > 0. For each u G OM there is a solution 
hi : G — > R to the equation 



For n > 1, f gOdg = f gdg where dg is the Haar measure normalised to 
be a probability measure and O any matrix in G. The integral of ge with 



dx\ = x e t gle dt, dx t = H(x € t )(gle )dt, 




For any smooth function F : OM — > R with compact support, 



F(xt) = F(uo)+ [ dF(xl)(H(xl)gle )ds. 



Ahi(u, g) = dF(u) (H(u)ei) (ge , e*). 
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respect to the Haar measure on G vanishes. For G = SO(n) it follows also 
from J SO(n) geodg = J s „-i sds, see Proposition 3.2.1 in Krantz-Parks IfTTH . 
Denoting by Dihi and D 2 hi the differential of hi with respect to the first 
and the second variable respectively, 



hi(x e t,g £ t ) 



l y r t l f 1 

--h(uo, I) + -j=J2 / ( D 2hi) { x % , gl) (g e s A k )dw k s + - / Ahi(x e s , g £ s )ds 
V e k Jo e Jo 

t 

(DihiXs*^ (H(u e s )g e s e ) da. 
Plug this back to F(xf) to see that 

F(x\) =F(u ) + eJ2 (hi(x\,9l) ~ h(u ,I)^ - V^Yl I \ D 2hi) m , g%) (A k g^dw k s 
i k,i ^° 

-eV / (Dxh^e^iHiuDg^ds. 
i J o 

The tightness of the law of {x\ } can be proved similar to that of the pre- 

e 

vious theorems, c.f. Lemma 13.21 Furthermore for s < t the following 
convergence holds in L 1 , 

~ e j 3 ^ Dlhi ^9l)^ H ^ ue s)9 e s^o)ds 

e 

^ ~ J s VdF(H(xl)ej,H(x e a )ei)ds J (ge^e^A^ige^eijd^g), 

where /j is the unique invariant measure for the A diffusion on G. The 
proof is similar to that of Lemma 13.31 taking into account of the following 
computation 

-(DxhiX^g) (H(u)(ge )) = -VdF (H(u)ge , H(u)ei) A' 1 (ge ,ei) 

= -^VdF (H(u)ej,H(u)ei) {geQ,e j )A~ x {geQ,e i ). 
j 

Now we assume that A = |A L . We may assume that {A k } is an or- 
thonormal basis of g and A = and let 



a i,j 



-J (geo^^i^A 1 ) 1 {ge ,e i )dg, 



where g is the Haar measure on G. For i ^ j the cross term a i: j van- 
ishes. There is an element O G G such that Oei = — and Oej = e 3 -. 
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Furthermore J2 t Af = —^1 and J2i9^f — _2 V"5'^- h is eas Y to see 
that (iA L )- 1 (^e , e<) = -^(#e , e*). The integral / (ge , etfdg is indepen- 
dent of i. In fact / (ge ,ei) 2 dg = j (ge ,Oei) 2 dg, for any O G G. Since 
/ Si (d e o, Ci) 2 dg = 1 it follows that 

A f 2 4 



n — 1 J G (n — l)n 

Finally we see that 

^2 ^ dF (H(u)ej, H(u)ei) (ge , e j )A' 1 {ge Q , e>) 



h3 



A n A 

> VdF(f) u (uei), f)u(wei)) = — A#. 

' « in — In 



(n — 1 m (n — l)n 

i=i 

The two operators A# and A are intertwined by n, and x\ converges to a 

€ _ 

Brownian motion. □ 
4.3 Another Intertwined Pair 

At this point we discuss a question asked to me by J. Norris. Since the 
process on the orthonormal frame bundle encodes the Riemannian metric 
we expect to see the Riemannian metric manifesting itself in some form, 
e.g. in the form of the corresponding Laplacian operator. Does the system 
below have a non-degenerate limit which is not necessarily associated to 
the given Riemannian metric on Ml In general the intertwined system 
would look like the following, 

du\ = CH{u\) o db\ + —=H(u £ t )V(x £ t , g £ )dt + -^=A* k (u £ t ) o dw\ + -A* (u £ t )dt 

v e v e e 

dx\ = Cut o db\ + --U^, gt)dt. 

Below we compute a simple case. The argument, with suitable adjustments, 
remains valid for the general case. 

Example 4.4 For simplicity consider R™ x SO{n) with the standard connec- 
tion, and the SDE 

dgl = -^=9 £ t A k o dw k t 

^ 1 (4.7) 

dx\ = 8g £ o db t + -^g £ V(x £ t , g £ )dt. 
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Here V is a R n valued function such that J G gV(x, g)dg = where dg is the 
Haar measure. For example take V(g) to be a function of even powers of 
g. We assume that V is suitably bounded with its partial derivatives in x 
suitably bounded. The parameter 8 is to be chosen. 

Letting A* k (g) = gA k . C = \ Y.k( A lf> assume that it is \A L . Taking 
5 = y/e, formal computation by multi scale analysis shows that : 

Claim. The limiting law for x\ is govern by the partial differential equa- 
tion on R n : 

= - / L gV(Xj9)dx Co 1 (gV(x, g)p) dg, (4.8) 

where the integral is with respect to the Haar measure on SO(ri). 

If 5 = 1 it ought to have, in addition, a Am term on the right hand side: 

dp f 

tjt = Am ft — j L u v( X ,u)£o 1 (L u v( X ,u)P)dv(u), 

which we do not discuss rigorously. A drift term in the g equation can also 
be added. Another interesting regime to consider is £\ big\ o db\ instead 
of g\ o db t with <5j takes values from {1, y/e}. In this case, a non-Laplacian 
like equation would follow. In the case that Si are all equal and V(x, g) 
is independent of x, the system can be interpret as an intertwined pair 
through time scaling. 

Equation Q4.8D can be deduced by the methodology below. Let / : M — > 
R be a smooth compactly supported function and A AI the Laplacian on M. 
Then 



f(x e t ) = f(x )+S j df(g e s db s ) + -S 2 j A M f(xl)ds + — j df(glV(xl,gl))ds. 



o 



If h is solution to £ /i(x, g) = df x (gV(x, g)), then 

-7= / df(g e s V(x e s , g e s ))ds = y/eh(x e t , u e t ) - \/eh(x , u Q ) - y/eS / d x h(x e s , g e s )g e s db s 
v e Jo Jo 

- [ d g h{glA k dw k s )-^~eS 2 [ A M h{x^gl)ds 
Jo Jo 



L gt v£dM x t,gl)ds. 

Since S = y/e, it is now easy to observe that {x e t } is a tight family. Since g\ is 
a fast ergodic motion and x\ does not move much as t — > 0, under suitable 
conditions, 

Efjxj) - f(x Q ) Efjxj) - f(x ) 
limhm = limhm = L n va hixn). 

e^0t->0 t t^Oe^O t 9 * 

has the required limit. 
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4.4 Perturbation to Horizontal Diffusions 

Let M be a compact connected n-drmensional smooth Riemannian mani- 
fold with connection V. Let zu be the corresponding connection 1-form on 
the orthonormal frame bundle OM with Lie group G, where G is taken to 
be 0(n) or SO(n) depending whether M is oriented. The horizontal bundle 
is integrable when and only when the curvature tensor of V vanishes. The 
Lie brackets of two fundamental horizontal vector fields will in general con- 
tribute to a vertical motion. However perturbation to horizontal flows can 
still be discussed and in this case we should consider not its projection to 
the manifold M unless the connection V is flat, but its motion transversal 
to the holonomy bundle. 

Let m G OM and r : [0, 1] -> M be a C 1 curve with r(0) = tt(m ). Let 
f be the horizontal lift of r through u Q . The parallel displacement fi of 
u can be written as u a some a e G. The set of such a that represents 
parallel displacements of u forms a subgroup of G and is called the holon- 
omy group with reference point u e OM which we denote by $(u ). In 
another word a e $(u ) if u and u a are connected by a horizontal curve. 
Denote by $o(wo) the restricted holonomy group which contains only loops 
that are ho mo topic to the identity loop. By Theorem 4.2, in Kobayashi- 
Nomizu IfToll . $(wo) 1S a closed subgroup and a sub-manifold of G with 
$0(^0) its identity component. Since M is connected all holonomy groups 
are isomorphic. 

Two points of OM are equivalent if they are connected by a C 1 horizon- 
tal curve. For each u in OM let P(u ) be the holonomy bundle through u , 
it consists of all u G OM such that u ~ u , i.e. u and u are connected by 
a horizontal curve. We may consider OM as disjoint union of sets of the 
form P(u). Let H = $(u ), which acts on P on the right, and P/H be the 
modulus space of P with respect to the equivalent relation. Then P/H is a 
smooth manifold and it can be identified with the associated bundle with 
fibre G/H and the equivalent relation: (uh^ 1 , h£) where £ denotes the coset 
corresponding to the identity. We identify (u, a£) with the orbit in P/H that 
contains ua. Denote by Ilx : P — > P/H the natural projection. The main 
task in the proof of the theorem below is to make sense of freezing the 
conserved Variable' and averaging out the 'fast 'variable'. 

Theorem 4.4 Let M be a connected and compact Riemannian manifold with 
a Riemannian connection V. Consider 

m 

dut = H(uf) o db t + H (u e t )dt + y/e Z k (u £ t ) o dw\ + eZ Q {u e t )dt, 

t! (4.9) 
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where H is a horizontal vector field and Z k are vertical vector fields. Then 
Hi(u\) converges in law, which is identified in &4.10D below. Furthermore 

define g\ by u\ = x\g\, where x\ = 7r(uf) and x\ its horizontal lift. The 
projection of g\ to the space of cosets, G/$(u ), converges weakly. 

Proof By the holonomy theorem of Ambrose-Singer [ 1 ] the Lie algebra of 
$(n ) is a subspace of q and is spanned by matrices of the form tt v (wi, w 2 ) 
where wi,w 2 are horizontal vectors at T P and v G P(uq). If u ~ v then 
= $(V). Let H = <P(wo), a manifold whose dimension is denoted n . 
We define a distribution S on OM: S = {T u (P(u)) : u G OM}. It is of 
constant rank, n + n . This distribution is differentiable and involutive 
and P(u) is the maximal integral manifold of S through u. Note that the 
holonomy bundles are translations of each other: P(u a) = P(u )a, a e G. 
If u is equivalent to v, the maximal integral manifolds through them are 
identical. 

Let u t be the solution starting from u of the equation 



Then u t is constant in t. To see this let / be a BC°° function on P/$(w ) 
and denote by III : P — > P/$(n ) the projection. Then 



/([«*]) = /(bo]) + / df (TUtiHius)) odb s + / df (TUMus)) ds. 



By the Reduction Theorem, page 83 of Kobavashi-Nomizu Ifl6ll . each holon- 
omy bundle P(u) is a reduced bundle with structure group and the 
connection in OM is reducible to a connection in P(u). Hence 



In particular we have TU^HTOM) = and f([u t ]) = f([u ]). 

We have shown that the solution to the horizontal SDE stays in P(u ) 
for all time. The horizontal SDE, restricted to the maximal integrable mani- 
fold P(u ), satisfies the Hormander conditions and is ergodic with a unique 
invariant measure Hp( Uo )- Fix a point u G P with x = vt(m ). Let v be the 
Haar measure on H = $(m ). Denote by v a the Haar measure on $(w a). 
Note that if v = u a some a G G, let u G $(w ) and a horizontal curve a 
with a = u , «i = u g, g G $(«o)- Then (3 = a a is horizontal with /3 = v 
and fix = u ga = va~ l ga. Consequently $(u a) = ad(a _1 )$(M ). If a G H, 
$(w a) = $(«o) and v a = v. 

Denote by N u the following fibre of the holonomy bundle P(u ): 



du t = H(u t ) o db t + H (u t )dt. 





T u (P(u)) = HT u OM © VT u (p(u)). 



N u = 7T~ (x) n P(U ), 



x = ir(u). 
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Locally N uoa = M x {ad(a _1 <E>(/uo)} and /xp (uoa) = dx x dv a where dx is the 
volume measure of the manifold M. 

Let F : OM ->■ R be a SC 1 function, we the integral 

JP(u) 

is defined to be a number depending on a transversal of P(u). On each 
fibre of the holonomy bundle P(u ) we choose a reference element v(x), 
which determines reference elements on holonomy bundles P(ua), due to 
that v(x)a is an element of P{ua) where u e 7r _1 0)- For any u e P(u ) there 
is g e H such that u = v(ir(u))g. We define 



Fd/i P{uo) := / F{v{x)g)dv(g)dx. 

P(uo) J M J P(u )r\ir- l (x) 

The resulting number is independent of the choice off. To see this let v' be 
another choice then v' = vh some h e H and u = va = v'h~ 1 a. Since G is a 
compact group, the Haar measure is bi-invariant, 

/ F{v{x)g)dv x {g) = [ Fiv'ix^^dv^g) 

F(v'(x)a')du x (g'). 



L 



' P(u )rm- 1 (x) 

Similarly if u = v(x)ag e P(u a) the following integral is well defined: 

/ Fdfi P(u) ■= F(v(x)ag)du a {g)dx. 

J P(u) J M J P(u)rnv-Hx) 

Evaluate / : P/$(w ) ->■ R at u\, denoting IIi(u) by [u], to see that 

/(K]) = /(K]) + V~e fdf (TU, (Z fc K))) o dw* + e /* d/ (TTIx (Z„K))) tfe. 

Jo Jo 

Let m be the Lie algebra of if and let Aj, i = 1, . . . , n be an o.n.b. of m. 

Let Bj, j — n + 1, . . . , N be an o.n.b. of the vertical part of the distribution 

S at u . Define A, = m U0 (Bj) e g. Consider the family of fundamental 

vertical vector fields j > n }, restricted to P(u ). Then TiIi(A*) = 

for i < n and for j > n , TJl^A*) = A*([u]). 

Writing Z k in terms of the basis {A k }, Z k = ^\ <J k A*, we have 

/(K])=/(K]) + V^^ E / ^K)d/W([<]))od«; 8 * 

j, — ii Jo 



fe i=no+l ' 



+ e V / ^«)d/(^(K]))tfe. 
„■_„. , i Jo 



j=no+l 
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The process [u e t ] is in general not Markov. It is however clear, follow- 
ing the standard method as used earlier, that the probability distributions 
{[wl],e > 0} is tight and any sequence of [u\] has a convergent sub- 

sequence with the same limit. The limit can be identified below. Define 



OijiM) = Efe>i / 0>fc^p( M) , and Z = Eln n +i a o A *j- For hj > n o define 



j=n +l u 0^j' 



^o(M) = / f °"o + o E da k^ Z ^> ) dftPM- 

J P(M) V 1 k>l J 

Then 

N 

£/([«])= ^ aij(N)Vd/(A;,A*) +df(Z([«])). (4.10) 

i,i=no+l 

The rest of the proof for the convergence is similar to that of Theorem 14.11 
For the second statement consider the process g\ defined by u\ = 7r(u e t )g^. 
Let p : G — > G/H be the canonical homomorphism. For a G G denote by 
[a] the left coset of H that contains a. Finally note that, if [u] denotes an 
element of of P/H that contains u, 

[u € t ] = Mu$)]gl = [u ]g £ t , 

and the motion in [u e t ] is essentially the motion of g\, while the latter is 
considered to be a representative of G/$(m ). Specifically, 

m N N 

dgl = V~tJ2Yl ^MdtiAj o dw\ + e viiulKAj, 

h=l j=l j=l 
m N N 

k=l j=n +l j=n +l 

The second statement thus follows. □ 

According to Theorem 8.2 in Kobavashi-Nomizu IfToTl if OM is connected 
there is a connection on OM such that P(u) = OM. On the other extreme if 
the curvature vanishes the orthonormal frame bundle OM foliates. This is 
so for a Lie group with the Left or right invariant connection. If M is simply 
connected the curvature zero case corresponds to the product bundle with 
the trivial connection. See also section 6.2 in Elworthy-LeJan-Li [11] for 
a discussion on the equivalence of the stochastic holonomy and holonomy 
and Arnaudon-Thalmaier [2] for work on Yang-Mills Fields and random 
holonomy. 
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